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ABSTRACT 
The  self-excitation  of  the  ballooning  tearing  mode  in  toroidal 
systems  with  strong  magnetic  fields  is  discovered.  This  instability  is 
driven  by  the  pressure  gradient  at  small  islands  and  arises  from  the 
coupling  of  poloidal  harmonics.  This  instability  is  always  unstable 
with  respect  to  perturbations  of  finite  size.  In  our  opinion,  it  is 
the  main  mechanism  for  plasma  turbulence  and  therefore  defines  the 
energy  confinement  time  in  tokamaks  and  stellarators. 
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In  every  contemporary  auxiliary  heating  experiment  in  tokamaks, 
the  loss  of  energy  is  anomalously  large  independent  of  the  level  and 
method  of  heating.  Although  the  invariance  of  the  loss  mechanism  is 
indicative  of  the  magnetohydrodynamic  nature  of  the  transport,  the 
physical  mechanism  for  this  loss  has,  until  now,  not  been  explained. 
We  will  show  that  in  a  high  temperature  toroidal  plasma,  in  which  the 
condition  B*Vp  is  well  satisfied,  a  finite  size  MHD  perturbation  of  the 
ballooning  tearing  mode  type  is  always  unstable  and  is  driven  by  the 
local  pressure  gradient  near  a  small  magnetic  island.  The  formation  of 
magnetic  islands  arising  from  the  ballooning  tearing  mode,  in  our 
opinion,  is  the  dominant  cause  of  anomalous  transport  in  toroidal 
systems. 

To  investigate  the  instability  of  the  tearing  mode,  we  examine  the 

1  p 
simple,  widely-used  model  for  tokamak  configurations  '  when  the  metric 

p  p 

tensor   is  given  by  g-,  i  =  1  ,  g-]  2  =  0,  g22  =  a  and  goo  =  (R-a  cos  6)  . 

2 

In  accordance  with  the  theory  of  nearest  equilibrum,   we  will  calculate 

the  saturated  state  of  the  instability  by  examining  the  stability  of 
the  perturbed  equilibrium  with  an  island  of  finite  width,  w.  The 
dominant  effect  of  the  perturbation  is  the  flattening  of  the  pressure 
profile  within  a  distance  w  of  the  resonance  surface.  This  is 
consistent  with  the  numerical  simulation  in  [4].  The  flattening  of  the 
pressure  results  from  the  equilibrium  condition  B'Vp  =  0  and  is  valid 
whenever  the  island  formation  and  growth  is  slow  in  comparison  to  the 
propagation  of  sound  waves. 

For  simplicity,  we  model  the  effect  of  the   islands  by  setting 

p  (a)  =  0  for   |a-a  I  _<  w.   Thus  our  actual  calculation  reduces  to  the 

I 

determination  of  the  stability  of  the  step  function  profile  of  p  shown 
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in  Fig.  1.  Saturation  corresponds  to  marginal  stability  of  the 
modified  profile.  Since  the  rotational  transform  is  approximately  zero 
in  the  island,  the  island  induces  a  spatial  variation  in  the  shear  as 
well.  Near  the  separatrix  the  flux  averaged  shear  decreases  from  a 
small  constant  value,  S^  to  order  S  W.  The  spatial  variation  of  the 

shear  S(x)  enters  only  in  higher  order. 

m(a-a^) 

We   introduce   the  dimensionless  parameter  x  =  where 

^m 
nq(ajjj)  =  m  and  examine  a  trial  function  E,   of  the  form  e 


ax,e.c)  =  ei^'"^-"^)  I   Fj^(x)ei^e  (1  ) 

2, 


The  standard  expression  for  the  potential  energy  of  a 
magnetohydrodynamic   plasma    is 

roo  2 

6W   =    I  dx[[(il-S(x)x)Fj,]'2  +    [i-S{x)x]^Fl   +    [^  +   aUJF^ 

where  S(x)  is  the  shear,  U  is  the  magnetic  well,  U  =  -  (l  -  ^-J  and 

R      q2 

a(x)  =  -8irp  Rq  /Bg  is  shown  in  Fig.   1. 

We  assume,  for  simplicity  that  S,a  <<  1  and  that  the  instability 
is  localized  near  a  single  resonance  surface,  a.^.  This  last  assumption 
implies  that  destabilizing  effects  of  order  a  exp  (-1/S)  will  be 
neglected.  We  restrict  to  the  case  when  the  probe  function  of  the 
energy  principle  consist  of  only  three  poloidal  harmonics,  Fq,F^,F_^. 
The  Euler  equations  for  (2)  are 


dx        dx  o        0^2       o 


a(P+D  ) 


»  ^   a  (x)D 


(3a) 


(3b) 


D   -  D  =  —^  +  "^^^  F^ 
4      2    o 


(3c) 


where  P  - and  D  - Equations  (3b)  and  (3c)  have 

solutions  given  by 


^  =  iTJ 


a(t)FQ(t)e 


-  x-t 


dt 


i^) 


D  =  1 


:a(t)A(x-t)  .^]F^(t)e-|x-ti 


where     A(x)  -  1     when     x  >   0     and     A(x)  -  -1    for  x  <  0.      After  applying 
these  expressions   for  P  and  D  the  variational  principle  becomes 


6W 


4   J 


S(x)2[(xFq)'2  +    (xfQ)2   +    aUF^] 

a'(x)FQ(x)dx   j"^   [a(t)A(x-t)    +  .^^il]FQ(t)e"  I  ^"^  I      (5) 


We  insert  the  simple  trial  function  Fq(x)  =  e  l^l/x  corresponding 
to  the  tearing  mode.  For  the  step  function  model  of  the  pressure 
distribution,  6W  simplifes  to 


6W  =  I    S(x)2  (xfo^'^  ■"  ^^^q]  *  2  J_  aUF^dx  +  a^FCw)  ch  w  [j"  FQ(t)e  ^  "  i  e  ''Fq(w)] 

~2w     2  — ?w 

<  2S2  .  ^"""    .  ^  [-.577  -  iln  2w  -  1/(4w)]  [l-Ll_)e-2-  (6) 

—         WW  J  -.    2    ' 


Clearly  for  any  values  of  the  parameters,  .577  +  iln  2w  +  l/(4w)  =  0 
determines  a  bound  on  the  width  of  the  island.  This  bound  has  the 
numerical  value  W  =  .3.  In  these  units,  the  distance  between  resonant 
surfaces  is  1/S,  Not  only  are  all  nested  flux  surface  equilibria 
unstable  to  finite  size  perturbation,  but  the  combination  of  many 
different  poloidal  and  toroidal  harmonics  will  result  in  island 
overlap.  The  existence  of  equilibria  with  a  finite  number  of  islands 
is  therefore  doubtful.  Thus  the  plasma  must  consist  of  large  regions 
of  ergodic  field  lines  or  constantly  move  between  a  statistical 
ensemble  of  equilibrium  states. 

In  the  limit  as  W  tends  to  zero,  the  energy  is  clearly  negative 
and  the  magnetic  well  canot  stabilize  this  instability.  The  width  of 
the  island  is  determined  by  the  marginal  stability  of  the  pressure 
distribution,  i.e.,  6W  =  0.  The  island  width  is  displayed  as  a 
function  of  a  for  fixed  U  and  S  in  Fig.  2.  When  a  is  much  less  than 
/8  S  the  island  is  small  and  satisfies 


W-  -^ (7) 


U  +  /U^  +  S^/2 


The  substitution  of  aj^^W/m  for  W  converts  this  expression  to  a 
dimensional  width. 
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Surprisingly,  the  quasilinear  calculation  has  shown  the  appearance 
of  a  new  linear  instability  on  the  modified  profile.  Furthermore  the 
limit  as  W  tends  to  zero  is  discontinuous.  Any  finite  size  island  in  a 
toroidal  system  will  cause  this  instability  even  when  the  plasma  is  far 
from  marginal  instability  of  linear  theory  for  nested  flux  surfaces. 
We  call  this  phenomenon  plasma  self -excitation. 

The  instability  is  caused  by  the  large  second  derivative  of  the 
pressure  near  the  separatrix,  p  (agg  ),  driving  the  satellite 
harmonics.  In  the  traditional  tearing  mode  in  a  cylinder  the  larger 
derivatives  of  j(a)  and  p  (a)  are  destabilizing  on  the  inside  of  the 
island  and  stabilizing  outside.  Thus  these  effects  cancel  to  leading 
order.  However  in  a  toroidal  geometry,  the  satellite  harmonics  are 
destabilized  on  both  sides  of  the  island.  This  effect  is  independent 
of  the  assumption  that  m,n  are  much  greater  than  one. 

The  calculation  reveals  the  presence  of  a  lower  energy  state  with 
islands  and  6W  is  the  release  of  potential  energy.  Thus  our  6W  can  be 
viewed  as  the  generalization  of  the  resistive  energy  principle  of 
Furth  to  a  toroidal  geometry  with  islands,   6W  can  also  be  interpreted 

r 

as  a  scaled  A  for  tearing  ballooning  modes  with  islands. 

Implicit  in  our  calculation  is  the  assumption  that  the  island  is 
larger  than  the  resistive  layer  and  thus  we  expect  the  island  to  evolve 
in  the  Rutherford  regime.  For  thermonuclear  plasma,  the  resistive 
layer  is  extremely  small  and  so  effectively  there  is  no  nonlinear 
threshold  for  the  size  of  the  perturbation.  Since  for  all  values  of 
beta,  there  are  equilibrum  states  with  islands  at  lower  potential 
energy,  the  assumption  of  nested  flux  surfaces  is  always  questionable. 
Furthermore  the  overlap  of  islands  occurs  at  much  lower  levels  of  beta 
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than  previously  believed.  Therefore  we  hypothesize  that  the 
destruction  of  flux  surfaces  resulting  from  the  tearing  ballooning  mode 
is  a  principle  cause  of  anomalous  transport. 

Our  stability  criteria,  AW  >  0,  determines  whether  the  potential 
energy  outside  the  resistive  layer  drives  the  instability.  Thus  our 
results  are  independent  of  the  model  of  dissipation  used  in  the 
resistive  layer. 
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APPENDIX:   TWISTING  MODE 

We  will  now  calculate  the  nonlinear  behavior  of  the  odd  twisting 
mode  in  physical  space.  We  exactly  reproduce  the  results  of  Zakharov, 
Semenov  and  Riedel,  however  our  present  calculation  is  substantially 
simpler  and  more  intuitive  than  the  previous  analysis  in  Fourier  space. 
We  again  base  our  analysis  on  the  theory  of  nearest  equilibrium  and 
therefore  analyze  the  stability  of  the  nested  flux  surface  equilibrium 
with  a  pressure  distribution  given  in  Fig.  2. 

Our  starting  point  is  the  energy  function  given  in  Eq.  (2).  The 
analysis  of  the  twisting  mode  requires  the  determination  of  the  next 
order  terms  and  the  retention  of  the  second  satellite  harmonic.  The 
general  Euler  equations  are 


^  ('n-Sx)2  i.  F,  -  (m-Sx)2F, 

=  f^  ^  [_  .  aU)F^  -  -  [F^,i*F^_i)  ^  -  [-^^  -  -^^]    ^  -   [ g j   (Al) 


This  equation  can  be  solved  with  a  Green's  function  to  yield 


e-|x-t|f 

F„(x)  =  -  -     -, —  (A2) 

">        2  J-a>  (m-Sx)(m-St) 


The   following  identities  will  be  useful 


f"    dx 


dF^    r^  e^'-^rjt)  [xS(x)]'f^(x) 


X  (m-Sx)(m-St)       (m-Sx) 


F„  -  

"1    dx 


dF^    rx    et-^f^(t)     [xS(x)]'f^(x)  ^^^^ 


-OD  (m-Sx)(m-St)      (m-Sx) 


We  assume  the  standard  ballooning  mode  ordering,  that  a  is  small, 

2 
the  shear  S  and  the  magnetic  well  W  are  order  a  and  that  the  width  of 

the  island  is  also  small.  We  neglect  the  destabilization  due  to  the 

resonance  of  the  satellite  harmonics  at  their  corresponding  resonant 

surfaces.  This  assumption  of  the  localization  of  the  eigenfunction 

around  the  resonance  surface  of  the  fundamental  harmonic  is  completely 

equivalent  to  the  averaging  of  the  equation  over  small  shear  in  the 

Connor-Hastie-Taylor  formalism.   The  inclusion  of  other  resonance 

surfaces  in  this  analysis  would  result  in  the  presence  of  the  standard 

destabilizing  term  of  order  a  exp  (-1/S)  in  Eq.   (A7). 

2 
With  these  assumptions,   F^-j   are  order  a  and  F^2  ^""^  order  a  • 

Thus  f^2  *^^"  ^®  approximated  by 


±2  -      2    *  dx    2     dx   4 


For  convenience  we  divide  F^^  into  the  sum  of  F°i   and  F^^   where 
F°  is  forced  by  the  fundamental  harmonic  F^.  Thus 
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n^ 


2  ~  dx   2  '^   dx  iJ 


2 

,1   _  "  ^±1    a  „   -  d  f°'^±2^   da  ^±2 
±^  -   2     2   ±2   dx  ^  2  ^    dx   4 


Integrating  Eq.   (A2)  by  parts  for  F°^  and  F°^  yields 


+  1  -  2  J 


a(t)F,et-x  .  ) 


^F^(t)e-|x-t| 
-00  dt  " 


-1  - 


P^.  -jJ__»<t,P,e-t-.  J_^^,^,,„-|x-t 


(A4) 


Furthermore,   we  can   integrate   (A3)   by  parts  to  find 


r  dF^pi        a   (x)F^,  „ 

a(x)    F^2   ±  ^     ~  ^^  *    0^«   W) 

±^         dx       -  8 


We  use  this  expression  in  (A3)  to  determine  the  desired  combinations  of 
^1 


dF 


,x-t 


«(>^)[^1  ^-n^]  = 


+li    a(x) 


dx 


[af,  - 


!^  .  £^  F„/2 
dt     dt  0 


7  2^  1 


(1-Sx)(1-St) 


+  a(xS(x)]'F^^(x) 


ct''(x)F 


o  ^  a(x) 


[[^  -    [tS(t))'s)F^  -  I  a^FAt)]e^-^   ^   a(xS(x)  ] 'f,.  (x ) 


■dt 


o   4  "  M 


-1 1 


dF 


a(x)[F_i  -  ——] 


It    a(x) 


dx 


^  o     dt     dt  0 


1  a^FAt)] 


(1+Sx)(l+St) 


+  a[xS(x))'F_i (x) 


"  ^""^^o       a(x) 


fx 


H       i 


[[^   +  2(xS(x)]']f„  -  1  a^F_At)]e^-''   +  a[xS(x) ) 'f_,  (x) 


■dt 


o   ^ 


These  expressions  can  be  simplifed  by  using  the  approximate  identity 


e^  ^aa^(t)F_, (t)  ^  ^ 

—  00  2   J 


e^"V(t)    e'^""a(u)F.(u)du 
J  X 


1 
2  J 


e^"^  a^  J   et-^d(u)Fo(u) 


'A. 


e^   ^a3(t)FQ(t)  +  a^Cx) 


e^  "a(u)FQ(u)] 


where  a  is  average  value  of  a(x). 

We  are  now  able  to  eliminate  the  satellite  harmonics   and 
re-express  f^  as 


^o  =  °'"f^o  "  -§- 


1  fda- 


£jF^(t)A(t-x)-|x-t|  .;  [Z^-  [xS(x))' 
dt  ^  ^      5id 


■|x-t| 


aF, 


f^  dx 


a(t)A(t-x)FQ)e' 


■lt-x|  _1 


1 
2  J 


-  t-z 


If„i 


(A5) 


where  A  is  defined  in  Eq.  (4).  If  the  width  of  the  island  is  set 
equal  to  zero  our  equation  is  the  Fourier  transform  of  Strauss 's 
equation.  Until  now  the  calculation  has  been  valid  for  an  arbitrary 
distribution  of  pressure  a(x)  and  shear  S(x)  in  the  island  region.  We 
now  restrict  to  the  step  function  distribution  for  a(x). 

The  contribution  of  the  island  on  the  potential  energy  is 


■12- 


•^"isl  =  a2F(w)w[y  -  -^—^ ]  (A6) 


where 


V  - 


e"^F„(t) 
w 


and 


e"  -  e"" 
W(w)  = — 


Since  Fq(w)  has  a  logarithmic  singularity  this  energy  will  be 
negative  and  thus  the  islands  have  a  destabilizing  effect  on  the 
twisting  mode  as  well.  Thus  twisting  modes  also  exhibit 
self-excitation,  however  the  destabilizing  term  scales  as  the  width  w 
while  the  corresponding  term  in  the  tearing  mode  scales  as  1/w,  Thus 
the  presence  of  islands  destabilize  the  twisting  mode  only  near  the 
marginal  stability  boundary.  We  determine  the  stability  criteria  for 
twisting  modes  by  evaluating  6W  with  the  trial  function  Fq^^^  ^  '*'o^'^^' 
the  modified  Bessel  function. 

This  trial  function  is  the  exact  solution  for  marginally  stable 
twisting  modes  when  the  magnetic  well,  U,  is  set  to  zero  and  there  are 
no  islands.  Evaluation  of  (A6)  gives 


6W  =  (1-a2)  ]l!£!  .  ^  .  a^wK  (w)[l  -  !^]  (A7) 


P    2       7   2 
where  X  ^  __  [l  -  )  and  A  =  1  corresponds  to  marginal  stability 
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when  U  =  0.   In  an  Initially  marginally  stable  equilibrium  the  islands 
will  grow  to  a  width 


,0205 


In   Strauss 's   model   the   region  of  linear  instability  is 

/2    7    a^       /2  ? 

1  -  -—  <  -L-  -   —  <  1  +__  and  the  center  of  this  region  is  S/a  =  7/16. 
M    1 6    S        4 
2 
The  value  a  /S  =  16/7  not  only  has  the  maximal  value  of  the  linear 

increment  but  also  the  largest  island  width 


w  =  .02  +  .015S 

Finally,  the  potential  of  the  magnetic  island  is  at  a  minimum  when 
w  =  .006  and  gives  a  nonlinear  stability  criteria 


6W  =  (1-a2)I!£!.^-  .001a2 


If  AW  is  positive,  the  potential  energy  outside  the  resistive 
layer  is  stabilizing.  Thus  when  the  twisting  mode  is  in  the  region  of 
AW  >  0  and  A  >  0,  the  mode  must  be  driven  inside  the  resistive  layer. 
Hence  an  instability  of  this  type  will  saturate  with  an  amplitude 
inside  the  resistive  layer.  Furthermore  the  stability  of  this  region 
depends  on  the  model  of  dissipation  chosen.  Finally,  when  the  twisting 
mode  has  AW  >  0  and  A  >  0,  the  mode  tends  to  different  values  of  b^  on 
opposite  sides  of  the  mode  rational  surface.  Thus  the  eigenfunction 
must  be  considerably  altered  by  resistivity  to  produce  this  mode. 

We  note  that,  due  to  the  crudeness  of  our  model  of  the  pressure 
distribution,  these  results  are  only  qualitatively  but  not 
quantitatively  correct. 
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FIGURE  3 
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